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Abstract 

00 

We study the long-time asymptotics of the doubly nonlinear diffusion equation p t — 
div(|Vp m | p - 2 Vp m ) in R n , in the range < m < and 1 < p < oo where the 

' mass of the solution is conserved, but the associated energy functional is not displacement 

convex. Using a linearisation of the equation, we prove an Z^-algebraic decay of the non- 
i-G ' negative solution to a Barenblatt-type solution, and we estimate its rate of convergence. We 

then derive the nonlinear stability of the solution by means of some comparison method be- 
tween the nonlinear equation and its linearisation. Our results cover the exponent interval 
< p < 2 ™+l where a rate of convergence towards self-similarity was still unknown for the 
p-Laplacian equation. 
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In this work, we consider the doubly nonlinear equation defined for any (t, x) € (0, oo) x M. n by 

(1.1) 



-£ = A p p:=div \Vp m \ p - 2 Vp m , (ieK",f>0) 



dp 
dt 

p(t = 0) = po > , [x G W 1 



with l<j><oo,0<m and n > 3. This class of equations contains the linear diffusion 
equation, (p = 2, m = 1), commonly known as the heat equation, dtp = Ap ; the nonlinear 
diffusion equation dtp = Ap m , known as the porous medium equation (p = 2, m > 1), or 
the fast diffusion equation (p = 2, m < 1), and the gradient-dependent diffusion equation, 
dtp = div(|V p\ p ~ 2 V p) := A p p, that is, the p-Laplacian equation, (p ^ 2, m = 1). When p ^ 2 
and m^l, Eq. is called the doubly nonlinear diffusion equation, due to the fact that its 
diffusion term depends non-linearly on both the unknown density p, and its gradient V/9. Such 
gradient-dependent diffusion equations appear in several models in non-Newtonian fluids [23], in 
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glaciology [201 EL an d i n turbulent flows in porous media [23]. For more details on these models, 
we refer to the recent monograph of Vazquez [31j . and the references therein. 

Assuming that the initial data is integrable, po G L 1 (M ra ), it is known that (|l.ip has a unique 
solution p G C ([0,oo),L 1 (W 1 )), with p(t) € C 1 - a (R n ) for some a G (0,1), see for instance 
[13 EH El [2S] . Moreover, starting with a non-negative initial data, pq > 0, it is known that the 
solution p{t) remains non-negative at all times. Furthermore for n > 3, there exists a critical 
exponent, 

n — p 



n(p — 1) 

such that if m > m c , then the mass of the solution is conserved, f Rn p{t) dx = J R „ po dx, while if 
m < m c , the solution vanishes in finite time, see |13U31] and the references therein. In particular, 
for the p-Laplacian equation, this corresponds to the critical p-exponent, 

2n 

Pc ■= — — r, 
n + 1 



above which the mass of the solution is conserved, while the solution disappears in finite time if 
p < p c . Therefore, up to renormalising the mass of po to unity, we can assume without loss of 
generality that, under the condition m > m c , the solution p(t) of (jl.ip is a density in R n , for all 
times t > 0. 

By similarity and scaling, it can be shown that, above the critical exponent m c , Eq. (|l.ip has 
a unique self-similar solution po,, whose initial value is the Dirac mass at the origin, that is, the 
fundamental solution of Eq. (jl.lH . In fact, among all the radially symmetric solutions of (11. ip . 
this solution is the most concentrated whose initial data have the same mass as po- It is called 
the Barenblatt solution [3J, and it is precisely: 



where 
and 

ud, (y) 



5 P := n{p — l)(m — m c ) > 0, 



a \ p ) p-1 



mp l , p 



1' 



with a and D* are uniquely determined by the mass conservation: ||ud, IIl^r™) = \\pD„ (t) \\L 1 (R n ) = 

IIpoIIl^k 71 )- 

When p = 2 and m > 1 — 2/n, the existence and uniqueness of the Barenblatt solution 
was proved by Friedmann and Kamin in |19| . Moreover, they showed that the solution p(t) of 
the Cauchy problem converges to PD„{t) w.r.t. the L 1 (M n )-norm, as t — > oo, with no rates. 
Rates of convergence were computed by Carrillo and Toscani [9 J if m > 1, independently by 
Del Pino and Dolbeault [16], and Otto [29] if m > 1 — 1/n. The rates found in this range were 
generically optimal. In the range 1 — 2/n < m < 1 — 1/n, there were studies of the linearised 
problem by Carrillo, Lederman, Markowich and Toscani [10], and Denzler and McCann [12] . 
These linearisations were useful to obtain rates of decay for the nonlinear fast diffusion equation 
by Carrillo and Vazquez [11] and later by McCann and Slepcev [27], and Kim and McCann |21j . 
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The decay rates obtained by using the linearisations are in general non optimal and is optimal 
in some sub-range, see [2T] , 

When p ^ 2 and m = 1, Kamin and Vazquez [22] proved existence and uniqueness of the 
Barenblatt solution pjj t for the p-Laplacian equation when p > p c , along with an L 1 -convergence 
of the solution p{t) of the Cauchy problem to pD t {t), with no rates. Their proof extends to the 
doubly nonlinear equation as long as m > m c , see |31j . Rates of convergence were computed by 
Del Pino and Dolbeault [17] when p c + l/(n + 1) < p < n for the p-Laplacian equation, but their 
rates are not optimal; see also a similar result for the doubly nonlinear equation in [18] . In [TJ [3] , 
Agueh generalises previous results by deriving optimal rates for the convergence of the solution of 
the Cauchy problem (jl.ip to pD*{t), for all m > m c + l/{n{p — 1)) = (n — p+ l)/(n(p — 1)) and 
p > 1. For instance, when p = 2, this condition coincides with the case m > 1 — 1/n, while for the 
p-Laplacian equation (p / 2, m = 1), it corresponds to p > p c + l/(n + 1) = (2n + l)/(n + 1), 
and therefore covers the range p > n left in [T7], but not the remaining exponent interval 
2n/{n + 1) < p < (2n + l)/(w + 1)- Similarly, for the doubly nonlinear diffusion equation, the 
rate of convergence remains unknown in the range 



m c < m < m c + 



1 n — p + 1 

n{p — 1) n{p — 1) 



(1.2) 



Indeed, the proof of [3] is based on optimal transportation inequalities, which follows from the 
displacement convexity [2U] of the energy functional associated with (jl.l|) . that is, H F (p) = 
J a „ F[p] dx, where 



F(x) 



p — 1 



rax ' 



-xlnx 



if m 



7 = m H if m ^ 



p — 1 
1 



I. 7(7-1) p-1 "" ' p- 1" 

This energy functional is displacement convex if and only if 7 > 1 — ^ or equivalent m > 
(n — p + 1)/ {nij) — 1)). This explains why the method of [3] does not extend to the interval (jl.2p . 

The goal of this work is then precisely to derive a rate of convergence w.r.t the L 1 (M ra )-norm, 
of the non-negative solution p of the Cauchy problem (jl.ip . to the Barenblatt solution pD»{t), as 
f — ► 00, provided that m belongs to the remaining exponent interval (|1.2p . that is, 



n — p n — p+1 

< m < 



nip — 1) nip — 1 

For convenience we rewrite the Cauchy problem (jl.ip as: 

( dp 



(1.3) 



-£ = div{pVc* [v(F'op)] }, (ier,f>o) 



(1.4) 



k p(t = 0) = po, (x€R n ), 
where c*(x) = \x\ p /p is the Legendre transform of the convex function 



/ \ l^l 9 1 1 

c{x) = , - + - 

q p q 



1. 
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By rescaling in time and space p as follows: 



where 

r = lnR(t), V = ^L, R(t) = (l + 6 p t)V s >, S p = (p - l)(nm + 1) + 1 - n, (1.6) 

it is easy to show that p solves (II. 4p if and only if u solves the rescaled convection-diffusion 
equation 

|H = divjuVc* [V (F'on)] + uy} (y £l",T> 0) ^ ^ 

_ u(t = o) = p (t/e» n ). 

Moreover, by conservation of mass there exists a unique such that the Barenblatt profile 
Ub» is the equilibrium solution of (|1.7p . Remark that in the considered range of exponents, 
m(p — 1) — 1 < and the Barenblatt profile is simply given by 

i 

u D , (y) = [D* + — -c(y)) 7-1 . (1.8) 

In fact, uo, is the unique density function of same mass as uq which satisfies on its support, 

V (F'ou D , +c) = 0. (1.9) 

The main result of our paper is the following: 

Theorem 1.1 (Rates of convergence) Let m be in the range (ll.3p and uq a density such that 
there exist positive constants Dq > D\ for which 

UD (x) < Po(x) = u (x) < u Dl (x) Viel". (HI) 

Consider u a solution to (11 -T[> with initial data uq, there exists a unique such that u(t) 
converges to the Barenblatt profile ud„ in L l {W n ). Moreover, there exist a time tq and two 
positive constants A and M = M(m,n,p,UQ,To) such that, for any time t > tq 

\\u(t) -u D ,\\ L i m < Me~K (1.10) 

As a consequence, for a time large enough the corresponding solution p(t) of (jl.ip converges to 
the Barenblatt solution /9D,(t), algebraically fast in the L 1 -norm, at the rate X/(25 p ): there exist 
a time to and a constant C = C(m, n,p, po,to) such that, for any time t > to 

\\p(t)-p D At)\\LH^)<Ct~ X /^\ (1.11) 

where 5 P = {p — l)(nm + 1) + 1 — n. 

The main tool is the following relative free energy with respect to the Barenblatt solution 
UD t defined by 



£[u\u Dt ] := / [Fou{y) - Foud^v) - F' ou Dt (y)(u(x) - u D ^(y))] 



dy (1.12) 
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for any given u G L+(R n ). Its derivative along the flow of (jl.4p is formally given by 

-~^;£[ u ( t )\ u d*} =1[u(t)\u d J 

where 

1[u(t)\u d J := / u(r,y) V (F'ou(r,y) + c(y)) ■ (Vc*oVF'ou(t, y) + y) dy . 

In this paper, we prove that the relative entropy decays exponentially fast in the form 

£[u(t)\u d .] < e-" T f[«o !«£».], (1.13) 

for some (3 > 0. This is obtained in two steps. First, we linearise (|1.7p at the equilibrium solution 
up, by using the linear perturbation u(r) = ud„ + eu(r), and we show that the linearised 
version of the relative energy converges to exponentially fast, as in [TO]. For that, we use 
the Hardy-Poincare inequality recently established by Blanchet, Bonforte, Dolbeault, Grillo and 
Vazquez in [5]. Next, following the strategy in [6], we try to compare the relative energy and 
the dissipation of the relative energy -that is, the Fisher information- for both linearised and 
nonlinear equations, to deduce the exponential decay (|1.13p for the nonlinear equation. The 
main differences with respect to [6] lie in the fact that a direct relation between the linearised 
and the nonlinear Fisher information is not clear due to the singular characters at the origin of 
the weights when 1 < p < 2. Therefore, we are forced to use a sort of regularised linearised Fisher 
information instead. Moreover, the control of the additional terms appearing in the regularised 
entropy dissipation of the linearised problem and in the relation between the entropy dissipations 
is more involved in our case. 

We note that, based on our computations (see Remark l3.4p . the Bakry-Emery approach used 
in [TO], which consists of differentiation twice the relative energy E[u(r)] to estimate the spectral 
gap at the eigenvalue 0, does not yield a positive result for our equation when 1 < p < 2, and 
thus, a similar procedure to [11] for the doubly nonlinear equations is not feasible. Moreover, the 
Hardy-Poincare inequality used here to establish the linear stability is actually valid on a larger 
interval, m, < m < m c + n ^_i^ , which includes our interval m c < m < m c + n ^_i^ , as < m c , 

where m* := + |rf- Therefore, our linearisation result extends naturally to the interval 
m* < m < m c where mass conservation for the nonlinear equation fails. In this range, one needs 
to carefully define the right class of initial data and a substitute of the Barenblatt solution, as 
done in [6] when p = 2. Here, we will not follow this path and we will restrict ourselves to the 
case m c < m < m c + n ^_i^ where mass is conserved to concentrate in the main new difficulties. 

The paper is organised as follows. In Section [21 we review and introduce the main estimates 
on the solutions needed in the rest of the work. In particular, we prove the convergence of the 
solution u(t) of (|1.7|) to the equilibrium solution v,£> t in C 1 (IR n ), as r — > oo, with no rate. Then 
in Section [3l we analyse a suitable linearised problem for which we apply an entropy-entropy 
dissipation argument based on Hardy-Poincare inequalities. Finally, Section H] is devoted to 
establish the exponential decay of u(t) to uq, by the comparison between linear and nonlinear 
relative entropy dissipations. 
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2 Convergence without rate 



Let us start by reviewing some well-known facts about the global unique weak solutions associated 
to in the range of exponents considered. They conserve mass for all times, i.e., 



p(t,x)dx = / po(x) dx ViG (0,oo). 
From now on, is the unique positive real such that 

/ po(x) dx = / U£) t (x) dx . 

Moreover, solutions of the Cauchy problem to enjoy a comparison principle and the L 1 - 

contraction property. Due to the change of variables (jl.5p . these properties hold for the solution 
u of the nonlinear Fokker-Planck equation (jl.7|) . Since in the rest of this paper we will only work 
with the scaled solutions of the nonlinear Fokker-Planck equation (jl.7p . from now on we will use 
t instead of r and x instead of y for the time and position variables respectively. The quotient 
function 

. . u(t,x) 
w{t,x) : = rT 

is solution to 

(jU) I 

d\v {w{t,x)u D ,{x)[Vc*oVF'o[w{t,x)u D ,{x)]-Vc*oVF'ou D ,{x)]} . 



dt u Dt (x) 
Define 



rrr ■ r U D Q { X ) / U Dl (x) 

W := mf — u —- < sup — := W\ . 

XGR" U D * \X) ,61" U Dt (X) 



A straightforward calculation gives 

with strict inequalities unless po = U D*- in terms of wq = uo/ud*, the "sandwich" assumption 
on the initial data (HI) of Theorem 1 1 . 1 1 can be rewritten as follows: there exist positive constants 
Do > D\ such that 

0<^o<^44<^o(^)<^44<^i<°o VxelT. (HI') 
u Dt \x) u Dt (x) 

Remark 2.1 Let us point out that the condition (HI) or (HI') in the case of the fast diffusion 
equation (p = 2) and in the corresponding range, 1 — 2/n < m < 1 — 1/n, is not restrictive. In 
fact, as a consequence of the Harnack inequalities proved in [7j, the hypothesis (HI) is satisfied 
by p{t) for any t > with an initial data uo £ L^(M n ). In the present case, a similar Harnack 
inequality, not available in the literature, would restrict the study of the asymptotic rates to this 
particular set of initial data. 



As a consequence of the regularity theory of degenerate parabolic equations [13], we can 
control uniformly C 1,Q -norms. 
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Lemma 2.2 (Uniform (^'"-estimates) Given a solution u G C([0,oo);L 1 (W 1 )) of (jl.7p with 
initial data uq satisfying (HI), then for any to G (0, oc) ; 

sup lie 1 -" (r^) < 00 • 
t>t 

Moreover, there exists C > such that for any x G M n \ B\ 

| Vu ;(x)|<C^ (2.1) 
I a; I 

Proof. Due to the comparison principle and the hypothesis (HI), the function u(t) is sandwiched 
between the two Barenblatt profiles for all times, i.e., 

ud < u(t) < u Dl t > 0, 

and thus is uniformly bounded in B2, the Euclidean ball of radius 2, uniformly in t > to > 0- 
Due to the regularity theory of degenerate parabolic equations |X3|, [HI \TE[ [25] , interior regularity 
estimates in the sense of u(t) G C l,a (Bi), for any < a < 1 hold uniformly in t > to > 0. 
Consider w = u/urj,, then w is also bounded in C ,a (Bi), for any < a < 1 uniformly in 
t > to > 0. To deal with large values of x we introduce, the rescaled function 

u x {t,x) := XP^-^u^Xx) 

which is also solution to (|l-7j) but the annulus B2\/B\ gets mapped into the annulus fii := B2/ B\. 
Note that all derivatives of the rescaled Barenblatt ^_d,/a are uniformly bounded from above and 
below since 

a R ( g(P-l) \ 

D u d*/\ — * CD 1 I I a; I ™(p-i)-i J uniformly in Qi as A 

for any multi-index f3 G N rf . As a consequence, we get that w\(t) = u\(t)/u Dt /x is a l so uniformly 
bounded from above and below in £li uniformly in A > 1 and t > to > 0. Again using the 
regularity theory of the degenerate parabolic equation, we deduce that the C 1,Q -norm of w\(t) in 
ill is also uniformly bounded for t > to and A > 1 by a constant C\. Going back in the A-scaling 
we find a constant independent of A > 1 such that 

|Vw(i,Ax)| Ci . . 

' 77 \ f < -f in (t 0> 00 (2.2) 
w{t,Xx) A 

in Qi, and thus, the C 1 -norm of w(t) in M. n /B\ is uniformly bounded. Similar scaling argument 
applies to the Holder semi-norms. As a consequence of f)2.2[) 



00 



1 , . 1 „ w(t,Xx) „ wit.Xx) 

Vw (t,Xx) < Ci— 4 — - < 2Ci \ , ; . 

A A|x| 

We thus obtain the desired result for any y = Xx G B2X \ B\, and any A > 0. □ 

From this we can obtain the following result regarding the evolution of the relative entropy 
to the stationary state. 
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Proposition 2.3 (Entropy /entropy production) Let u G C([0,oo); L 1 (M n )) be a solution 
of (|1.7p for an initial data satisfying (HI), and consider the free energy £ defined by (|1.12p . Its 
derivative along the flow of (|1.7|) is: 

—£[u(t)\u D „] = -l[u(t)\u Dt ] 

where 

l[u(t)\u Dt ] := [ u(t) [V (F'ou(t) - F'ou Dt )] • [Vc* (V (F'o U (t))) - Vc* (V (F'ou Dt ))]dy 

is the relative Fisher information of u(t) w.r.t. «o t . Moreover, Z [u(t)\uD,] = if and only if 
u = U D, ■ 

Proof. By performing formally integration by parts, we get 

-^-£[u(t)\u Dt ] = [ [F'o U (t) - F'ou D *]div\u(t)Vc* [V (F'o U (t))]+u(t)y\ dy 

dt J R n I J 

= - / «(t)V [F'ou(t) - foujj.] • {Vc* [V (F'on(t))] + y\ dy. 

The above energy dissipation follows using that ud„ satisfies (|1.9p and Vc* o Vc = id. This 
integration by parts can be justified using Lemma 12.21 by a standard argument introducing a 
cut-off function like in [BJ Proposition 2.6]. Since the arguments are exactly equal, we do not 
perform any further details. By the convexity of c*, 

[Vc*(a)-Vc*(b)}-(a-b)>0 (2.3) 

with equality if and only if a = b. So the Fisher information is non-negative and zero if and only 
if u and ud„ have the same mass and such that V (F'ou(t) — F'ouaJ = 0, i.e. u = Ud»- D 

With these ingredients, we can obtain a first result of convergence toward stationary states. 

Lemma 2.4 (Uniform convergence) Let u G C([0, oo); L 1 (W 1 )) be a solution of (|1.7p for an 

initial data satisfying (HI), then lim^ 00 w(t,x) = 1 uniformly in compact sets ofM. n . 

Proof. Define u (t, x) := u(h + 1, x), for any given h > and t € [0, 1]. It is also well-known [13] 
that equi-bounded set of solutions of (jl.ip are equi-continuous in time. This property carries 
over to u(t) by the change of variables in (jl.5p . This fact together with the uniform bounds in 
C 1,a obtained in Lemma [2 .21 and the Ascoli-Arzela theorem implies that for any sequence (/i n )neN 
there exists a sub-sequence (/i n )neN 5 denoted with the same index, such that {u hn } ne N converges 
to a function Uoq uniformly in compact sets of [0,1] x M n , and moreover, u 00(f) € C 1,a (R Tl ) for 
all t G [0, 1]. Since £[u(t)\uD t ] is non-increasing and positive and 

rh„ + l rl 

£\u(h n )\u Dx ] -£[u(h n + 1))|«dJ = / l[u(s)\u Dt ] ds = / 2[u(s + h n )\u Dt ] ds , 

Jh n JO 

the function t 1— > I[u hn (t) \ud„] is integrable on [0, 1] and converges to zero as n — > 00. By (12. 3p . 
X is non-negative. By Fatou's lemma we have 



lim u hn (t,x) V (F'ou hn (t,x) + c(x) ) ■ (Vc* VF'ou hn (t,x) + Vc*oVc(x 



dx = . 
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As a consequence of (|2.3[) . Mqo satisfies V[-F'ouoo(x) + c(x)] = 0, from which Uoo = ud for some 
D > 0. By conservation of mass D = D*. Since the limit of all the convergent sub-sequences is 
uniquely determined by ud 9 , the result is proved. □ 

Proposition 2.5 (Convergence in L p -spaces) Let u € C([0, oo); L 1 (]R n )) &e a solution of the 
scaled doubly-nonlinear equation (|1.7j) /or an initial data satisfying (HI), i/ien lim^oo ||tt(i) — 
«£»,||p = 0, for any p G [l,oo]. 

Proof. By Lemma [2.41 lim t ^oo \u(t,x) — ud,(x)\ = for any x £ M n . Moreover, by assump- 
tions (HI), for |x| large enough 

\u(t) — ud, I < max { \ud — ud* \ , |""Di — "p, [ } = O ^|x|~9(2-7)/(i-7) 

So the difference between \u(t) — UD t f is in L 1 (M n ) if 9 > Q(p,m) with 

n(l - 7) 



&(p, m) 



/(2-7) 



It is easy to check that Q(p,m) is a decreasing function of 7 and so of m. Since q > 2, in the 
range of exponents (jl.3p . we have 

0(p, m) < B(p, m c ) = < 1 . 

n + q 

By Lebesgue's dominated convergence theorem, it implies that u(t) converges to ue>„ in L (M. n ), 
for any 6 £ [1, 00). Finally, we use the following interpolation lemma, due to Nirenberg, |28} p. 
126]: 

n 2 

||/||oo<C||/||^ ( 2 Rn) ||/||r 2 V/eC(R"), 
for / = u(t) — ud* together with Lemma 12.21 to obtain the result in the uniform norm. □ 

Remark 2.6 In contrast with [6J, we do not generally have the convergence in C 1,a (W ri ). 



3 Linear stability 

To prove the decay (|1.13p of u(t) to ud„ in the energy form, it is sufficient to establish the 
following logarithmic Sobolev type inequality: 

£ [u\u D *] < -I[u\u D *] , (3.1) 
P 

for some (5 > and u € C 1,a (IR n ) verifying (HI). Indeed, (|3,ip combined with Proposition 12.31 
yield 

[u(t)\u D *] < -PS [u(t)\u D *], 

at 

and this leads to the energy decay (j!.13[) by a simple Gronwall argument. To prove (|3.ip . we will 
first show a linearised version of this inequality, by considering the perturbation u(t) = ud* +ev(t) 
of the solution u(t) to (II. 7)1 . This will lead to the convergence of v(t) to in relative entropy for 
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the linearised equation of (jl.7p . as we will show below. Next section will be devoted to compare 
the relative entropy and Fisher information in (|3.ip with their linearised analogues. 

For clarity in our exposition, let us start by formally deriving the form of the linearised 
logarithmic Sobolev inequality that we will be dealing with below. Using the perturbation u = 
up + ev and the second order Taylor expansion of F (ud* + ev) at e = on the expression of 
the relative entropy (|1.12j) . we have that 

Fou - Fou D * = evF'(u D *) + —v 2 F"ou D * + 0(e 3 ), 
and then £ [u\uo*\ linearises as: 

£ [u\u D *] = — / v 2 F"o UD * + <3(e 3 ). 
2 it" 

Let us hence introduce the weighted L 2 -norm: 

E[v} = ~[ v 2 (x)F"ou D *(x) dx, (3.2) 

which will play the role of the linearised relative entropy. 

Concerning the linearisation of the Fisher information, from the first order Taylor expansion 
of F'ou = F' (ud* + ev) at e = 0, we have 

B ■= V [F'ou] =A + eW + 0{e 2 ). (3.3) 

with 

A := V [F'oud*] = Vc and W := V [vF" ou D *] . 
Then using that Vc*(z) = z\z\ p ~ 2 , we obtain 

Vc*(.B) = S/c*(A) + e\A\P- 2 W + e(p - 2)\A\ P ^{A ■ W)A + 0{e 2 ). (3.4) 

Combining (13. 3p and (13. 4p . we see that Z[u|ii£i*] formally linearises as: 

X [u\u D »] = e 2 [ u D * \A\ p - 2 \W\ 2 dx + e 2 (p-2) [ u D * \A\ p - 4 (A ■ W) 2 dx + (9(e 3 ). 

Hence, for e small enough, the logarithmic Sobolev inequality (|3.ip linearises as 

v 2 F"ou D *<( u D *\Vc\ p - 2 \V[vF"o UD *]\ 2 dx (3.5) 

+ (p-2) / u D *\Vc\ p - A [Vc-V (vF"ou D *)] 2 dx. 

It will be shown below that the l.h.s of (|3,5p is a Lyapunov function - and the relative entropy 
- for the linearised equation of (jl.7p . and the r.h.s of (|3.5p corresponds to the dissipation of this 
relative entropy, up to a constant. 

Let u be the solution of (jl.7p . and consider the small perturbation 

u(t) = u D * + ev(t) (3.6) 
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of «d*, where e > is small, and v(t) G C 1,a (W l ) for some a G (0,1). Because of the mass- 
conservation, we have that 



/ v(t,x)dx = 0, Vt>0. 



Moreover (|3.6p implies that 

du dv 

dt = e m- (3 - 7) 

On the other hand, using (jl.9p . we have that Vc* (A) = Vc*[— Vc(x)] = — x, and then (|3.4p gives 
that 

u [Vc* (JB) + x] = enzj. [|yl| p - 2 Ty + (p - 2)|^| P " 4 (A • + 0(e 2 ). (3.8) 

Inserting f|3.7[) - ()3.8[) into (jl.7p . we formally obtain after simplifying by e and then setting e = 0, 
that the linearised problem to (II. 7p is: 



^ = divlu D * (\A\p- 2 W + (p - 2)U| p - 4 (A • WOA) 1 (xel",i> 0) 

a* i j (3.9) 

(x G R n ), 



with wo G L 1 (M n ) of zero average. We can easily check that Eq. (|3.9p has the linearised relative 
entropy (|3,2p as Lyapunov functional. Actually, differentiating E [u(t)] along a solution v of (|3.9p . 
we formally have by a straightforward computation, that 

-^E [v(t)] = - (I [v(t)] + (p- 2)l [«(*)]) , (3.10) 

where 

\[ v (t)]= [ \W(t)\ 2 u D *\A\ p - 2 dx and \ [v(t)]=[ (A ■ W{t)f u D * \A\ P ~ A dx. 

The Cauchy-Schwarz inequality implies that |^4| p (A ■ W(t)) < \A\ p ~ 2 \W(t)\ 2 , and as a conse- 
quence, lo [v(t)} < I [v(t)]. Using 1 < p < 2, we have 

I [v(t)] + {p- 2)l„ [«(*)] > (P- 1)1 [«(*)] > 0. (3.11) 

In case p > 2, it is direct that I [t> (i)] + (p — 2) lo [u(<)] > I [v(t)] > 0. From these estimates, the 
dissipation (|3,10p and |A(x)| = |V [F'oud*] (x)| > for all x G M. n , we readily formally conclude 
that the unique steady state is the zero solution. 

The objective of the rest of this section is to show the following asymptotic exponential 
relaxation of the linearised equation (13. 9p : 



Theorem 3.1 Let m satisfying (|1.3p and vq G L 1 (M. n ) with zero average. Consider v(t) the 
solution to (|3.9p with initial data vq. There exists a constant (3 > such that 

E[v(t)]<e-^E[v ]. (3.12) 

Let us concentrate first in the case 1 < p < 2. To derive this exponential rate of convergence, 
we establish the following functional inequality 

E[v]<^(\[v] + (p-2)\ [v]), /3>0, 
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for all v E C 1,a (M. n ) with zero average. This inequality corresponds to the formal linearisa- 
tion (|3.5p of the logarithmic Sobolev inequality (|3.1|) . In fact, because of (|3.1ip . it is sufficient to 
prove the following linearised logarithmic Sobolev type inequality: 

EH<^IH, /?>0 (3.13) 

for all v £ C 1,a (W l ) with zero average. This is equivalent to show the Hardy-Poincare type 
inequality: 

/ g 2 dfi(x) < f \Vg\ 2 dv(x), (3.14) 

if™ JR™ 

for some > 0, where g is any function satisfying f Rn g dfj,(x) = 0, and 

d/i(x) = - dX du(x) = u D *(x)\V (F'o UD *( x )) r 2 dx. (3.15) 

Note that the functions v and g in (|3.13p - (|3.14j) are related by g = v F"oud*, and (3 = 2(p — l)//3. 
The inequality (|3. 14[) is also enough to prove the needed inequality in the case p > 2 since f|3. 10|) 
implies that 

■^E[ W (t)] <-![«(*)]. 

Lemma 3.2 (Hardy-Poincare type inequality) Let m,n,p be such that 1 < p < oo and 
m c < m < • Then, there exits a constant (3 > such that 

g 2 d/i(x) < p I \Vg\ 2 du(x), 

for any function g G C 1,a (IR n ) satisfying J^„ g d//(x) = with < a < 1, where fi and v are 
defined by (|3.15p . 

We keep calling this inequality, "Hardy-Poincare inequality" to remind the link with the 
inequality proved in [3 E] but here we only use the Poincare type part of the inequality. The 
proof of the Hardy-Poincare inequality was performed in [5] and an estimate of the constant ft is 
also established. This proof can be adapted to Lemma 13.21 For completeness, we give here the 
proof of this variant of the Hardy-Poincare inequality. 

Proof. We first observe that we can reduce to show the inequality for the Schwartz class g £ D(R n ) 
by simple approximation arguments. In L 2 (M n ,d/i) we consider the closable quadratic form 
v h- > Q(v) := L n |Vg| 2 du and — £ the unique non-negative, self-adjoint operator in L 2 (M n ,d//) 
associated with the closure of Q. By Persson's theorem {301 



inf o~ css (— C) = lim inf 



_R— >oo vgHr f Rn \g\ 2 dfi 

where TLr := {v G H 1 (R d , dv) : supp(w) C W 1 \ B(0,R)}. Roughly speaking it means that 
the inequality is true for any weights with the same behaviour in a neighbourhood of +oo. By a 
straightforward computation using (jl.9|) and (|1.8|) . we have that 

dfi(x) = -(d* + 1—1\x\A 7 " dx -i^i^oo - (1—1) \x\ 2a ~ 2 dx 



12 



and 

du(x) = \x\ 2 - q ( D* + — ^\x\ q J dx ~i a .i_ +00 ^) |x| 2a dx 

\ mg / 11 \ mq J 

with a chosen in such a way that g(2 — 7)/(7 — 1) = 2(a — 1) and 2 — q + q/i'y — 1) = 2a, that is, 

5(2-7) , , 2-g « 
« = 1 + 777 7T or equivalently a = — h 



2(7-1) " 2 2(7-1) 

It is left to the reader to check that a < — (d — 2)/2 in the range of 1 < p < 00 and m c < m, and 
thus, we can apply Theorem 1] to obtain 

info- ess (-£) > - — -«a 

where k q is the constant of the following Hardy inequality with weight. We refer to [5j Theorem 
1] for estimates on the constant K a depending if — n < a or a < —n. We remark that both cases 
happen depending on the precise values of m c < m < n^p-i) an< ^ ^ < P < 00. 

The lowest eigenvalue of — C is Ai = with eigenfunctions given by the constants functions. 
Zero mean- value solutions belong to the orthogonal set to the eigenspace associated to Ai. Since 
Ai is non-degenerate we obtain the desired result for some A2 £ (0, K a ]. □ 

Since the only behaviour of the weights that counts is their growth at infinity, we can avoid 
the singularity of the weight at the origin for the singular case 1 < p < 2 directly to obtain the 
following stronger inequality. 

Corollary 3.3 For any e > 0, there exists a constant f3 e > such that 

g 2 dfi(x) < k f \Vg\ 2 du e (x), 

for any function g € C 1,a (IR n ) satisfying L n g dfi(x) = 0, where [i is defined by (|3.15p . and 

du e (x) = u D * (x) (e + I V [F'ou D * (x)) \) p ~ 2 dx. 

Therefore, setting v = g/(F" oud*), we have the stronger linearised logarithmic Sobolev in- 
equality 

E[v] < |/« [v] (3.16) 
for all v E C 1,a (R n ) with zero average and < a < 1, where 

[v}= [ \V(vF"ou D *) | 2 (e+ |V (F'ou D *(x)) \) p ~ 2 u D * dx. (3.17) 



Remark 3.4 The Bakry-Emery approach used in |10j to establish the linearised logarithmic 
Sobolev inequality (13.13h when p = 2, does not seem to apply here when 1 < p < 2. For illustra- 
tion, let us consider the particular case m = n = 1, that is the linearisation of the 1-dimensional 
rescaled p-Laplacian equation, dtv = div{(p— 1)ud* \A\ P ~ 2 W}. In this case, the relative entropy 
dissipation equation ()3.10p simplifies as 



dt 
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and it is easy to show that its dissipation is 



dl[v(t)] 



2(p-lfV[v(t)\v D ,] 



dt 



where 



V[v{t)\v D *] = 



F"o UD * [di-v (u D *\A\ p ~ 2 



W)f dx. 



Following [10], if one can establish the estimate T> [v (t)\vD*] > \l[v(t)\vD*] for some A > 0, 
then it will imply the linearised logarithmic Sobolev inequality (|3.13[) . But by a direct computation, 
we can show that 



Note that the second term in the above expression is non-negative (because F is convex), while the 
last term is non-positive in the range 1 < p < 2. If p > 2, the last term is also non-negative and 
we obtain T> (v(t)\vD*) > l{v{t)\vD*), that is A = 1. This then yields the desired inequality f)3. 13j) 
when p > 2. But if 1 < p < 2, one cannot derive this estimate from (13.180 . at least at cursory 
glance. 

Proof of Theorem^J\ We apply Lemma[3J2]to g(t) = v(t)F"ou D * and = 2(p-l)//3. By (IXTH 
and (I3.11|) we have 



4 Nonlinear stability 

The first step to go from linear to nonlinear stability is to use that our solution is sandwiched 
between two Barenblatt profiles to compare the nonlinear relative entropy and its dissipation 
with their linearised counterparts. 

Proposition 4.1 (Comparison linear /nonlinear relative entropy) Consider a function u 
satisfying (HI). Then there exist positive constants C\ and C% such that 




(3.18) 




This leads to (|3.12j) by a Gronwall estimate. 



□ 



C\ E[u — udJ < £N"UdJ < C2 E[u — ud,\ ■ 



Proof. By Taylor's formula on the integrand of the relative entropy we have 




with 



ud*(x)Wq < min(u(x),UD t (x)) < £(x) < max(«(x), uo t (x))) < ud,(x)Wi 
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due to (HI), see also (HI'). The asserted result follows from homogeneity of F" with C\ := 
mW 7 ~ 2 and C 2 := mW^' 2 . □ 

The next objective is to compare the nonlinear Fisher information, I \u{t)\uD*\, with its linear 
analogue, I \u{t) — ud*] along solutions of (II. 7|) . Let us point out that the weight 

\V(F'ou D ,)r 2 = \x\ 2 -i 

in the linearised entropy dissipation diverges at the origin for 1 < p < 2. This singular behaviour 
makes complicated any attempt to compare it with nonlinear Fisher analogues. Due to the 
singularity of the weight |V (F'oud*) [ p ~ 2 at x = 0, we will replace I \u{t) — up*] by its regularised 
analogue 1^ [u(t) — ud*] defined by (|3.17p . where e > is a fixed constant. 

Proposition 4.2 (Comparison linear/non-linear Fisher information) Assume that u is 
the solution of ()1.7j) . and set v = u — ud* . Then: 

1. Case 1 < p < 2: Given e > 0, there exist to > and positive constants k\ and K2 such that 
for all t > to, 

l (e) [v(t)} < Kl l[u(t)\u D *) + K 2 E[v(t)] . (4.1) 

2. Case p > 2: There exist to > and positive constants K\ and k 2 such that for all t > to, 

l[v(t)\ < KlI[u(t)\UD*]+ K 2 E[v(t)\ . (4.2) 

Moreover k 2 can be chosen arbitrary small provided that to is large enough. 

The proof of this proposition is organised as follows: 

Claim 1: We first show that, for all e > 0, there exists Ko > such that 

\U[v(t)]<K \U[v(t)}+K 2 E[v(t)], (4.3) 

where 



/ |V [F'ou- F'ou D *]\ 2 (e+ \V (F'ou D *(x)) \) p 2 u D * dx. (4.4) 
Claim 2: Next we show that if 1 < p < 2, then for all e > there exists 5 > such that 

\M[v(t)]<51[u(t)\u D *], (4.5) 
whereas if 2 < p < oo, then there exists S > such that 

l 7 [«(*)] < 61[u(t)\u D *) ■ 
Combining (|4.3p and (|4.5p . we obtain the desired inequalities (|4.ip - (|4.2p with k\ = 5 no- 
Proof of Claim 1: Here we follow the arguments of the proof of Lemma 5.1 in [6]. Indeed, let 
hk(w) = (w k ~ 1 — l)/(k — 1), where 

. . u(t,x) 
w(t,x)- 



u D *{x) 
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Because of assumption (HI), we have that Wq < w(t,x) < W\, where the constant Wo and W\ 
are such that < Wo < 1 < W\. By studying the function h^/h^ on [Wq, Wi], we have 



and 
where 

and 

Now, define 



ao/i 7 (w) 2 < h,2(w) 2 < a\h^(w) 2 , 



h' 2 (w) < a2h'(w) 



(4.6) 



a := I7 - 1| 



Wo-1 



W 



7-1 



1 



< 1, «i := I7 — 1| 



Wi - 1 



w 



7-1 



1 



> 1 



a 2 := W 2M > 1. 



JR™ V 



e + x 



q-l\p-2 



ud* dx. 



We have that I ^ = 1^ and for k = 7, l-y is defined in (|4.4I) . Next we compute \ k [v]. By 



expanding 



V (u' D ^h k {w) 



we have 



m l u 



h' k (w) 2 \Vw\ 2 (e + \x\ q - 1 ) p - 2 dx 



+(1 - 7 ) 2 / h k {w) 2 \x\ 2{ - q - 1 \e + [x^f^UD, dx 



'i e) H 



+2m(l- 7 ) / u' y D *ti k (w)h k (w)\x\ q - 2 (e+ \x\ q ~ l ) p ~ 2 Vw ■ x dx. 
Integrating by parts, the last integral can be rewritten as 
v? D *ti k {w)h k (w)\x\ q - 2 (e + {x^-y^Vw ■ x dx 

= 1 [ V (h k (w) 2 ) ■ \x\ q ~ 2 x(e + \x\ q - 1 )' p - 2 ul, dx 

2 ./ran 



\ f h k {w) 2 dw(\x\ q - 2 x(e+ Ix^-y^ul*) dx 
h k {w) 2 div (\x\ q - 2 x(e + {x^y- 2 ) dx 



+ 



2 
1 
2 

2m 



h^wflxl^-^ie+lx^-y^UD* dx. 



Then, 



\®[v] = m 2 I u 2 ^ 1 ti k {w) 2 \Vw\ 2 {e+\x\ q - l ) p - 2 dx 

+ (l- 7 ) f h k {w) 2 \x\ 2(q - 1) {e + \x\ q - 1 ) p - 2 u D * dx 

JR™ 

- m(l - 7) / /i fc (w) 2 div (|x| g - 2 x(e + \x\ q - 1 ) p - 2 ) u£. dx. 



(4.7) 
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Next we set kq := max(ai , a2 ) • Moreover, since w uniformly converges to 1 as t goes to oo, then 
ao, ai, a 2 and kq > 1 can be chosen arbitrary close to 1 provided that t > to, for some to large 
enough. Combining (|4.6|) - (|4.7|) . we have 

I (*)[«] = \®[ v ] <m 2 a 2 [ u%~ 1 h'Jw) 2 \Vw\ 2 (e+\x\ q - 1 ) p - 2 dx 

+ ai(l- 7 ) f h y {w) 2 \x\ 2(q - 1 \e+ ^^f^uo* dx 

JR™ 

- m(l - 7 ) f h 2 {w) 2 d\w (\x\ q - 2 x(e + |x| 9_1 ) p ~ 2 ) v? D , dx 
< koI ®[v] + m(l - 7 ) f u^div (|x| 9 - 2 x(e + |x| 9_1 ) p_2 ) (K h 7 (w) 2 - h 2 {wf) dx. 

Finally, we observe that < Koh^(w) 2 — h 2 (w) 2 < (ko/ q o ~~ 1) h 2 (w) 2 and by direct computation 

| div (|x|"- 2 x(e + [xl"- 1 )^ 2 ) | < n + 2(q - 2), (4.8) 

and we then deduce that 

I W [«(*)] <K \M[v(t)}+K 2 E[v(t)} 

with k 2 := 2(no/ao — 1)(1 — 7)(n + 2(g — 2)) > 0. Clearly k 2 is arbitrary small provided that to 
is large enough. □ 

Proof of Claim 2; case 1 < p < 2: First we expand l-y (i)] and Z [u(i)|u.D*], and we have that 
I 00 [«(*)] = / I V [F' o u] | 2 (e + | V [F' o Uz) .] |) p " 2 u D . dx 

+ / | V [F' o u D *] | 2 (e + | V [F f o u D *} \ f~ 2 u D * dx 

-2 [ V[F'ou] -V [F' o u d *] (e + |V [F' o u D »] |) p ~ 2 u D « dx 

and 

T[u{i)\u D *]= f \V [F' o u ]\ p u dx + f |V [F' o u D »] | p u dx 

7R" JR™ 

-/ V [F'on] • Vc* (V [F'ou D »])udx 

-/ V \F' o ud*] ■ Vc* (V \F' ou])ii dx. 
Jr™ 

Next we use Young inequality a ■ b < c(a) + c*(b) with c(z) = l-zl 9 /^, a = Vc* (V [F' o u]) and 
6 = V [F' o n£)»], to have that 

V[F'o UD «] - Vc* (VF'(n)) < i|V[F'o-u]| p + i|V[F'on D .] | p . 
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Then X [u(t)\u£>*] can be estimated as 

l[u(t)\u D *)>- f \S7 [F' o u ]\ p u dx + - [ |(V [F' o UD *])\ p udx 

P JR" '/ JR n 

-[ V[F'o«] ■Vc*(V[F'ou L) *])udx 

JR n 

= - I \V[F' o U ]\ p udx + - [ | (V [F'oud*]) \ p u dx 



nm 



dx. (4.9) 



7(1-7) 

Now, we compute the cross term of I7 We have that 

f V[F'o«])-V[F'o u D *} (e+| V [F' o u d *] \) p ~ 2 u d * dx 
Jm." 

V(u 7_1 ) ■ x\x\ q - 2 (e + \x\ q - l ) p - 2 u D * dx 



m 



1-7 



1 


- 7 




m 


1 


- 7 




1 


1 


-7 



/ — div (x\x\ q - 2 (e + Ix^-y- 2 ) dx 
JR" w x 

+ — !— [ w^- 1 u D ,\x\ 2( - q -^(e+\x\ q - 1 ) p - 2 dx. 
1 - 7 Jr" 

Since the last term in the above sum is non-negative, then using (I4.8p and the fact that 1 < p < 2, 
we can estimate ^[^(i)] as 

\\ e) [v(t)] < [ \W[F'ou]\ 2 (e + \V[F'o UD ,} \y~ 2 u D , dx 

jR n 

+ I |V [F' o u D *\ | 2 (e + I V [F' o u d *] |) p "\ d . dx 

jR n 

+ _2m_ f ^ div , |a; |,_ 2(€+ , |,_i )P _ 2 x dx 

1 - 7 y Rn W 

< / |V [f o U ]| 2 ( e + |V [F' oii fl ,] |) p -\ D . dx 
+ / |V [F' o u D *\ \ 2 (e + I V [F' o UD *\ \) p - 2 u D * dx 

JR n 



1.7 



| 2m(n + 2(q-2)) _ ]( 

1-7 Ui 



< / |V [F' ou]\ p § e (u,u D *) 2 ~ p u D * dx+ [ \V[F' o UD *]\ p u D * dx 

Jl™ ~ jR n 



where 



2m(n + 2(q-2)) f if 
1-7 ii« w 



, , . IV \F' o u 



e + |V[F'onz). 
18 



Using Lemma 12.21 and (|2.ip . the reader can check that $ e (u,UD*) is uniformly bounded by 
expressing the gradients in terms of derivatives of w and ury*, and thus, & e (u,ur)*) < f] for some 
7] > depending on e. Therefore, we obtain the estimate 

\i e) Ht)] < tj- I \V \F' o u]\ p u dx + f \V \F' o u D *]\ p u dx 
+ 2m(n + 2(g -2)) y 



(1 - 7 )W 

Combining (|4.9|) and (|4.10p . and setting 

A 1 / 2-p 2 7 (n + 2(g-2)) 

«5:=— max pr/ 

Wo V n 

we deduce (|4"3j) . □ 
Proof of Claim 2; case 2 < p < oo: As above, we have the expression 

l 7 [u(t)]= / |V [F , o- u -F , on D »]| 2 |V(F'o-u D *(x)) | p " 2 n D . dx. 
it" 

For convenience, we can also rewrite I (u(t)\uj)*) as 

T[«(i)|« c *] = / H[u\u D *) |V [f on — F'oud 
Jk 



\ 2 u dy, 



where 

„. V{F'o U -F'ou D *)-[Vc*[^{F'ou)}-Vc*[V(F'o UD ,)]} 

H[u{t)u D * = — — — — 

| V [r'ou — p'oud*) \ 

and /C := {x € W 1 such that |V [F'ou — F'ouo*] \ ^ 0}. Let us show that there exist a constant 
5 > 0, such that for all t > to, 

H[u(t)\u D *] >5\V [F'oud*) \ p - 2 . (4.11) 

Let us remark, if p = 2, then 8 = 1, and equality holds in (|4.11|) . 

For simplicity, set a(t) = V(F'ou(t)) and ap* = V (F' our)*). It is clear that (|4.1ip holds 
in the set where dry* =0. Therefore, let us restrict to the set where ap* ^ without loss of 
generality. Let us denote b(t) = a(t)/\aD* \ and br>* = od*/|od*|. It is straightforward to check 
that 

H[u(t)\u D *} _ (b(t) -b D *)- (\b(t)r 2 b(t) - bp*) 
\V (F'ou D *) \p- 2 \b(t)-b D *\ 2 ' 1 ' ' 

Let 9 denote the angle between b(t) and bp*. We have that 

\b - b D *\ 2 = \b\ 2 + \b D *\ 2 - 2bcos9 = 1 + \b\ 2 - 2bcos6, 

and 

(b-b D *) ■ {\b\ p ~ 2 b - b D *)) = \b\ p - \b\ cos 9 - |6| p_1 cos 9 + \bp* | 2 = 1 + \b\ p - {\b\ + \b\ p - 1 )cos9 
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so that (|4.12p reads as: 

H[u(t)\u D .] _ l + r(t)P - (r(t) + rjtyP- 1 ) x(t) 
\V(F'ou D *) \P~ 2 ~ 1 + r(t) 2 - 2r(t)x(t) 

where r(t) = \b(t)\ > and x(t) = cos# G [—1, 1], with r(t) — > 1 as t — ► oo. Estimate (|4,lip is 
reduced to show that 

for all x G [—1, 1] and for all r > 0. For that, let us define the function 

F p (r, x) =:= 1 + r p - (r + r p ~ l )x - 5(1 + r 2 - 2rx), 

which is easily checked to be decreasing in x for r > 0, whenever 5 < |. Therefore, we have 
F p (r, x) > F p (r, 1) and thus, to show (|4.13p for r > and — 1 < x < 1 is reduced to show that 
Fp(r, 1) > 0, whenever 5 < ^. Since F p (l, 1) = 0, this is equivalent to show that f p (r,l) > 5 for 
< r < 1 and r > 1. The last assertion comes from the fact that when p > 2, 

(r p - r) - (rP- 1 - 1) _ (r - l)^ 1 - 1) _ r p ~ l - 1 
M r ' 1) - ( r - 1)2 - ( r - 1)2 ~ r- 1 

is bounded below by 1, since in < r < 1, we have r p ~ l — 1 < r — 1 < 0; in r > 1, we have 
r p_1 — l>r — 1>0; and lim r ^i 1) = p — 1 > 1. Therefore, 

1 

2 

□ 



/ p (r, !)>!>-> 8. 



Proof of the main theorem, Theorem Given 1 < p < 2 and e > 0, set v = u — ud*- From 
Proposition 14.21 and the strong linearised logarithmic Sobolev inequality (|3.16p . we have that 

EKt)]<-^iKt)M, 

2 - K 2 Pe 

where k 2 can be chosen arbitrary small provided that t > to is large enough. This together with 
proposition 14.11 yields the logarithmic Sobolev type inequality: 

£[u{t)\u D *]<^I[u{t)\u D *], (4.14) 
where A := ^2 — K2/3 e ^ /C2^i(3e > 0. We combine ()4.14p and the entropy dissipation equation 

-^S[u(t)\u D *} = -l[u{t)\u D ,} 

to obtain the exponential decay of the relative entropy, £ [u(t)\uD*] < e~ xt £ [uq\ud*\. The L 1 - 
decay ([l.lOp follows from the Csiszar-Kullback type inequality (see for e.g., [3]), 

\\u[t) - wd*|Ili(]r«) - M(n,n,p)£ [u(t)\itD*] , M(m,n,p) > 0, 

and (jl.lip is a direct consequence of the rescaling (|1.5j) - (|1.6|) . The case p > 2 follows analogously 
without need of using the regularised entropy dissipation. □ 
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Remark 4.3 The rate of convergence A in Theorem li.il can be explicitely reconstructed in the 
above computation for a given choice of e for p < 2, but to will depend on the choice of e. For 
p > 2, we can even give an explicit range for the constant A. 
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